Time-calibrated phylogeny of the woody Australian genus Hakea (Proteaceae) supports multiple origins of insect-pollination among bird-pollinated ancestors†
The authors thank A. B. Thistle for editorial advice, M. Simmons and two anonymous reviewers for comments on the manuscript, W. Maddison, M. Holder, and P. Midford for conversations about the methods, and individuals and institutions who provided field and laboratory assistance, leaf material, unpublished data, images, and primer suggestions, including J. Carter, D. Feller, G. Jordan, A. Kubes, R. Makinson, N. Miles, P. Olde, and G. Sankowsky. This research was funded by an NSF grant (DEB-0516340) to A.R.M. and by support from the Royal Botanic Gardens and Domain Trust, Sydney, to P.H.W. and from the South Australian Department of Environment and Natural Resources to R.M.B. and W.R.B.
Abstract
• Premise of the study: A past study based on morphological data alone showed that the means by which plants of the Australian genus Hakea reduce florivory is related to the evolution of bird pollination. For example, bird pollination was shown to have arisen only in insect-pollinated lineages that already produced greater amounts of floral cyanide, a feature that reduces florivory. We examine a central conclusion of that study, and a common assumption in the literature, that bird pollination arose in insect-pollinated lineages, rather than the reverse.
• Methods: We combined morphological and DNA data to infer the phylogeny and age of the Australian genus Hakea, using 9.2 kilobases of plastid and nuclear DNA and 46 morphological characters from a taxonomically even sampling of 55 of the 149 species.
• Key results: Hakea is rooted confidently in a position that has not been suggested before. The phylogeny implies that bird pollination is primitive in Hakea and that multiple shifts to insect pollination have occurred. The unexpectedly young age of Hakea (a crown age of ca. 10 Ma) makes it coincident with its primary bird pollinators (honeyeaters) throughout its history.
• Conclusions: Our study demonstrates that Hakea is an exception to the more commonly described shift from insect to bird pollination. However, we note that only one previous phylogenetic study involved Australian plants and their honeyeater pollinators and that our finding might prove to be more common on that continent.
Hakea Schrad. & J.C.Wendl. (149 species; Proteaceae) is a genus of shrubs and trees endemic to Australia. It reaches its greatest diversity in the fire-prone communities of Australia's southwest corner, where 92 species are endemic (5, 6). Species also are established now outside their natural range as invasives in similar environments in South Africa (37), New Zealand (79), and southern Australia (6). Hakea is strikingly diverse in ecologically significant traits, including flower form (specialized for insect, bird, and mammal pollination), leaf form (terete to broad, flat leaves), strategies for persistence in fire-prone communities (both vegetative resprouting and serotinous fruits), and strategies for nutrient and moisture acquisition below ground, and it has therefore proved useful to comparative ecology (see, e.g., 29; 6; 61; 30). Nevertheless, until now we have been without a well-supported phylogenetic framework for understanding the evolution of its diversity. Here, we report a phylogeny for Hakea inferred from molecular and morphological data, estimate the genus’ age, and explore the phylogeny's relevance to previously published conclusions about the direction of evolutionary shifts between pollinator classes and the evolutionary interplay between pollination and florivory.
Conventional wisdom in pollination biology holds that specialized bird-pollinated taxa have repeatedly evolved from generalized insect-pollinated ancestors (Faegri and van der Pijl, 1966; 23; 28; 15; 60). Moreover, evidence for the reverse transition from bird to insect pollination is rare. 73, for example, list 29 phylogenetic studies of pollinator shifts and associated morphological changes, of which 13 document a total of more than 35 transitions from insect to bird pollination. In contrast, only two of these studies (59; 77) found shifts from bird to insect pollination (two transitions from hummingbird to bee pollination in Sinningieae and five transitions from hummingbird to hawk moth pollination in Aquilegia). The classic evolutionary explanation for such an unbalanced trend is that specialization may be almost irreversible because it restricts “evolutionary potential,” though this idea is difficult to test empirically (24). The birds featured in all but one of these studies belong to the New World family Trochilidae, the hummingbirds, even though at least seven other bird families have been implicated in floral pollination (62). Thus our knowledge of the evolution of bird pollination is strongly biased to the highly specialized hummingbirds, and it is reasonable to ask whether the generalizations that have been made about the evolution of bird pollination really pertain only to them.
We know of only two phylogenetic studies of bird-pollinated systems from outside the New World that are comparable to those listed by 73, both involving the southwestern Pacific family Meliphagidae, the honeyeaters (13; 30). 13 analysis of the Nemcia-Brachysema-Jansonia clade (Fabaceae tribe Mirbelieae) showed two shifts from bee to honeyeater pollination, using a morphology-based tree as its phylogenetic framework. However, a subsequent molecular systematic analysis by Chandler et al. (2001) implied additional origins of bird pollination in this group, as well as the possibility of reversal to bee pollination in one clade. Secondary transitions from honeyeater to insect pollination seem more likely than analogous shifts from hummingbird to insect pollination because the degree of pollinator specialization found in many hummingbird-pollinated plants is not found in honeyeater-pollinated taxa (23). Most of the regular nectar-feeding species of honeyeaters visit a wide range of ornithophilous plant species, and most honeyeater-pollinated plant species are visited by a wide range of co-occurring birds. 30 found multiple origins of honeyeater pollination from insect-pollinated ancestors, with four to five reversals, in their study of Hakea (Proteaceae). However, like Crisp, they also used a weakly supported, morphology-based tree, that of 6, as the phylogenetic framework for their analysis.
In addition to reconstructing the history of pollinator shifts between insect and bird pollination, 30 examined the relationship between plant features relevant to pollination (e.g., stigma–nectary distance) and those relevant to discouraging florivory (e.g., floral cyanide concentrations and floral accessibility) in 51 species of Hakea in southwestern Australia. Damage to flowers through florivory and destructive foraging has the potential to greatly reduce plant reproductive success, and selective pressure for adaptations to deter such losses are likely to be high (51; 15), especially in resource-limited environments like those typically inhabited by Hakea species. This is especially so for plants pollinated by perching birds such as honeyeaters, the smallest of which are much larger, heavier, and stronger than the largest insect pollinators. 30 found evolutionary correlations between the mode of pollination (insect or bird) and the degree to which a species’ flowers are protected structurally and chemically (by cyanogenic compounds). Furthermore, they found evidence that “bird pollination could only have evolved from insect-pollinated taxa with a pre-existing capacity to protect inflorescences through the synthesis of floral cyanide” (p. 257). This finding is significant because theirs is one of only a few studies to show plant defenses “pre-adapting” lineages for pollinator shifts (76). 30 invoked a 50-Ma history to support their results, citing the occurrence of florivorous ancestors of the emu, cockatoo, and marsupials in Australia before diversification of the pollinating ancestors of the honeyeaters, a process that commenced within the last 45 Ma, after the origin of the superfamily Meliphagoidea in the mid-Eocene (3; 18) and before the late Miocene (8). 30 based their conclusions on analyses using a phylogeny inferred from morphological characters by 6 and the informal taxonomic groupings established by 5 for their Flora of Australia treatment. It is noteworthy that 6 consider bird pollination to be the more probable primitive condition of Hakea on the basis of the greater geographic sampling represented in their phylogeny (i.e., one not restricted in scope to southwestern Australia). Here, we revisit these conflicting conclusions, using a sampling of Hakea that is geographically wide and taxonomically even to produce a dated phylogeny inferred from morphological and molecular data.
Hakea was placed in a new subtribe, Hakeinae Endl. (tribe Embothrieae), by 75 with Grevillea R.Br. ex Knight (357 species), Finschia Warb. (three species), Buckinghamia F.Muell. (two species), and Opisthiolepis L.S.Sm. (one species). Grevillea is the third largest genus of angiosperms in Australia (55). Its monophyly has yet to be established (46), but that question is outside the scope of our study. Hakea appears monophyletic on the basis of the synapomorphy of woody follicles that dehisce down both sides of the style base; follicles are nonwoody and only split along one side in close relatives (34; 6). For their Flora of Australia treatment of Hakea, 5 established 31 informal groupings of species. These groupings were established before 6 inference of a phylogeny for the genus on the basis of 59 morphological characters from 88 species. 5 also were working with a strict 100-word limit on the morphological descriptions of groups, which led to groupings of more morphologically similar taxa than might have been otherwise circumscribed had it been easier to adequately address suspected homoplasies. 6 showed that some of the informal groups of 5 were not monophyletic, though not all groups were tested at the time. 6 named six clades comprising more than one informal group on their phylogeny—the largest, the Multilineata clade, was composed of six informal groups.
We tested the following hypotheses: (1) the clades of 6 are monophyletic, as are the informal groups of 5 not already tested by 6, and thus the groupings provide a useful framework for evolutionary studies such as that of 30; (2) Hakea arose before the Eocene, at a time when florivores were present, but the primary avian pollinators of Hakea (the honeyeaters) had yet to begin diversifying extensively, as implied by 30; and (3) bird pollination arose repeatedly from insect pollination in Hakea, as concluded by 30. We also considered the relevance of the rooting of 6 phylogeny to 30 conclusions about the evolutionary relationships between the pollen vectors used and floral cyanide concentrations.
MATERIALS AND METHODS
Taxonomic sampling
We created a set of molecular data consisting of four regions from the DNA of the chloroplasts and other plastids (cpDNA) and three nuclear DNA (nDNA) regions from up to 75 species. Data for 70 of these species are presented here for the first time; data from the rest first were presented by 48 or 49, 50). The 75 species include 55 (of 149) species of Hakea, one species from each of the other 11 genera in tribe Embothrieae, one species from each of the other three tribes in subfamily Grevilleoideae, one species from each of the other four subfamilies in Proteaceae, and one species from each of the other two families of Proteales (Appendix 1; following 75). Our sampling of Hakea was taxonomically even, drawing 1–6 representatives from each of the 31 informal groups recognized by 5. 30 assigned species to bird- or insect-pollination classes on the basis of stigma–nectary distances (SND; >13 mm for the former, <13 mm for the latter). According to this diagnostic, we conclude that we included both insect- and bird-pollinated representatives from the Ulicina group of 5 but only the more common pollinator classes from three other groups with species that fall into both classes (Lorea, Varia, and Multilineata; one species in each of these falls into the alternate class and was not sampled).
DNA extraction, amplification, cloning, and sequencing
We extracted genomic DNA using the DNeasy Plant Mini Kit (Qiagen, Valencia, California, USA). The four cpDNA regions we sampled were the matK, atpB, and ndhF genes and the rpl16 intron; the three nDNA regions were waxy loci 1 and 2 and the PHYA gene. We amplified and sequenced all the DNA regions except the rpl16 intron, following the procedures used by 50. Following 50, we cloned all the nDNA regions but not the cpDNA regions. We sequenced the rpl16 intron following procedures described by 49.
Defining substitution characters
We aligned and edited the DNA sequences using Sequencher software (Sequencher version 4.7 sequence analysis software, Gene Codes Corporation, Ann Arbor, MI USA). We compared the plastid sequences to the complete plastid sequence of Nicotiana tabacum L. (GenBank accession number NC 001879) to determine the beginning and ending point of each cpDNA region. We compared the waxy sequences to the complete waxy sequence for Solanum tuberosum (GenBank accession number X58453) to determine the beginning and ending points of introns and exons and the numbering scheme for the exons. Alignment of all the coding regions was straightforward and did not imply any changes in reading frames, only a few insertions or deletions that were multiples of three nucleotides in length. We aligned the rpl16 and waxy intron sequences using the criterion of similarity, without regard to optimizing some description of the inferred phylogeny (66). We excluded from the data rpl16 and waxy intron sequences from outside of subtribe Hakeinae because we considered them to be too divergent across order Proteales to be confident of the homology of aligned positions. For modeling of nucleotide substitutions in the data, we recognized each of the plastid regions, the PHYA gene, the waxy locus 1 exons, the waxy locus 1 introns, the waxy locus 2 exons, and the waxy locus 2 introns as separate data partitions.
Coding morphological characters
We (R.M.B. and W.R.B.) modified the original morphological data of 59 characters for Hakea (6), adding taxa not previously included. For the present analysis, we simplified two transformation series from those data: one for leaf shape and division was simplified from four characters to two (characters 5 and 6, Table 1), and one for pollen presenter shape was simplified from two characters to one (character 29, Table 1). Analyses with and without use of these transformation series produced very similar results. Of the remaining 56 characters, 46 are parsimony informative and were used in our study (table 2).
| Number | Character | Character states |
|---|---|---|
| 1 | Lignotuber | Absent (0), present (1) |
| 2 | Ferruginous hairs on branchlets | Present (0), absent (1) |
| 3 | Leaf base | Not twisted (0), twisted (1) |
| 4 | Leaf insertion | Not clasping stem (0), clasping stem (1) |
| 5 | Adult leaf shape in cross section | Flat (0), terete (1) |
| 6 | Adult leaf division | Simple (0), compound (1) |
| 7 | Abaxial surface of terete leaves | Grooved (0), not grooved (1) |
| 8 | Adaxial surface of terete leaves | Grooved (0), not grooved (1) |
| 9 | Mucro numbers of flat leaves | One (0), several to many (1) |
| 10 | Margin of flat leaves with single mucro | Entire (0), crenulate (1) |
| 11 | Prominent marginal veins on leaves | Absent (0), present (1) |
| 12 | Palisade layer of leaves | Continuous around leaf and not broken by vascular bundles (0), broken by vascular bundles capped by lignified parenchyma (1) |
| 13 | Inflorescence position | Axillary (0), axillary and terminal (1), terminal (2) |
| 14 | Involucre around incipient inflorescences | Present (0), absent (1) |
| 15 | Developed shoot on inflorescence rachis at time of flowering | Absent (0), present (1) |
| 16 | Rachis | Simple (0), branched (1) |
| 17 | Bract subtending flower pair | Present (0), absent (1) |
| 18 | Number of flowers in inflorescence | Less than 10 (0), 10–30 (1), over 30 (2) |
| 19 | Flower color | White, cream, or yellow (0); pink to red, white and pink, or red-flushed (1) |
| 20 | Rachis length | Obscure or knoblike (0), terete and less than 15 mm (1), elongate and 15 mm or more (2) |
| 21 | Flower buds | Straight (0), curved (1) |
| 22 | Degree of fusion of perianth parts | Four free parts (0), two or more parts at least partly fused (1) |
| 23 | Perianth length | Below 3 mm (0), 3–6 mm (1), over 6 mm (2) |
| 24 | Relative size and shape of tepals | Equal (0), unequal (1) |
| 25 | Ferruginous hairs on perianth | Absent (0), present (1) |
| 26 | Nonferruginous hairs on perianth | Present (0), absent (1) |
| 27 | Pollen color | Yellow to white (0), pink (1) |
| 28 | Nectary gland thickness | Thin flap (0), thick strip (1) |
| 29 | Pollen presenter shape | Oblique (0), erect (1), lateral (2) |
| 30 | Pollen presenter cone | Covering whole of base (0), surrounded by basal flange or disc (1), absent (2) |
| 31 | Retention of seed | Retained on plant in closed fruit >1 yr (0), retained on plant <1 yr (1) |
| 32 | White wood zone in fruit | None (0), present or reduced (1) |
| 33 | Extent of pale wood zone in fruit | Extensive (0), not extensive (1) |
| 34 | Fruit axis at base | Straight to slightly oblique (0), strongly oblique or curved (1) |
| 35 | Fruit axis at apex | Straight to slightly oblique (0), strongly oblique or curved (1) |
| 36 | Spines on fruit | Absent (0), present (1) |
| 37 | Pusticules on fruit | Absent (0), present (1) |
| 38 | Corky outgrowths on fruit | Absent (0), present (1) |
| 39 | Horns on fruit | Present (0), absent (1) |
| 40 | Apiculum at fruit apex | Present (0), absent (1) |
| 41 | Position of fruit apiculum | Apical (0), lateral (1) |
| 42 | Shape of fruit valves | Straight (0), recurved (1) |
| 43 | Extent of seed cavity on valve face | Much less (0), almost filling (1) |
| 44 | Position of seed wing | Apical (0), down one side only (1), down both sides in full or part (2), encircling (3) |
| 45 | Distal ridge on seed | Present (0), absent (1) |
| 46 | Surface ornamentation on seed body | Smooth-ridged (0), with tall projections (1) |
| Character | ||||||||||||||||||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Species | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 | 31 | 32 | 33 | 34 | 35 | 36 | 37 | 38 | 39 | 40 | 41 | 42 | 43 | 44 | 45 | 46 |
| Hakea arborescens | 1 | 1 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 2 | 0 | 1 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 |
| H. archaeoides | 1 | 0 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | ? | ? | 0 | 0 | 2 | 1 | 2 | 1 | 0 | 1 | 0 | 0 | 0 | ? | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | ? | 0 | 0 | 0 | 1 | ? | 0 |
| H. auriculata | 1 | 1 | 0 | 1 | 0 | 0 | ? | ? | 1 | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 2 | 1 | 0 |
| H. baxteri | 0 | 0 | 0 | 0 | 0 | 0 | ? | ? | 1 | ? | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 2 | 0 | 1 | 1 | 0 | 1 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 2 | 1 | 1 |
| H. brachyptera | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | ? | 0 | 0 | 0 | 1 | 0 | ? | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | ? | 1 | 0 | 0 | 0 | 1 | 3 | 1 | 1 |
| H. bucculenta | 0 | 0 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 2 | 1 | 2 | 1 | 0 | 2 | 0 | 0 | 1 | 0 | ? | 1 | ? | 0 | 1 | ? | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | ? | 0 |
| H. clavata | 1 | 0 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 2 | 0 | 1 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 2 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 2 | 1 | 0 |
| H. commutata | 1 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 |
| H. conchifolia | 1 | 0 | 0 | 1 | 0 | 0 | ? | ? | 0 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 1 | ? | 0 | 1 | 0 | ? | 0 | 0 | 1 | 0 | 0 | 1 | ? | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | ? | 0 |
| H. constablei | 0 | ? | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | ? | 0 | ? | 0 | 0 | 1 | 1 | 0 | 1 | 0 | ? | ? | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 2 | ? | 0 |
| H. corymbosa | 0 | 1 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 1 | 1 | 2 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 2 | 0 | 0 |
| H. cristata | 1 | 0 | 0 | 0 | 0 | 0 | ? | ? | 1 | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 2 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | ? | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 3 | 1 | 0 |
| H. cucullata | 0 | 0 | 0 | 1 | 0 | 0 | ? | ? | 0 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 2 | 1 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | ? | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 1 | 0 |
| H. dactyloides | 1 | ? | ? | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | ? | 0 | ? | 2 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | ? | 0 | 1 | 0 | 0 | 0 | 0 | 2 | ? | 0 |
| H. drupacea | 0 | ? | 0 | 0 | 1 | * | 1 | 0 | ? | ? | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 2 | 0 | 2 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | ? | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 |
| H. eriantha | * | ? | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 2 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 0 |
| H. grammatophylla | 1 | ? | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | ? | 2 | 1 | 2 | 1 | 0 | 2 | 0 | 0 | 1 | 0 | ? | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | ? | ? | 1 | 0 | 0 | 0 | 0 | 1 | ? | 0 |
| H. hastata | 0 | 0 | ? | ? | 0 | 0 | ? | ? | 0 | 1 | ? | 1 | 0 | 0 | 0 | 0 | ? | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | ? | 0 | ? | ? | ? | 0 | 0 | ? | 0 | 1 | 0 | 0 | 0 | 0 | 1 | ? | 0 |
| H. horrida | 0 | 1 | 0 | 0 | 1 | 1 | ? | ? | 1 | ? | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | ? | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 3 | 1 | 0 |
| H. incrassata | 1 | 1 | 1 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | ? | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | ? | 0 | 2 | 0 | 1 | 0 | 0 | ? | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 3 | 1 | 1 |
| H. invaginata | 0 | ? | 0 | 0 | 1 | 0 | 0 | 0 | ? | ? | ? | 1 | 0 | 0 | ? | 0 | ? | 1 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 2 | ? | 0 |
| H. lasianthoides | 0 | 0 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 2 | 2 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 0 | 0 |
| H. laurina | 0 | ? | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | ? | 1 | 0 | 0 | ? | 0 | ? | 2 | 1 | 0 | 1 | 0 | 2 | 0 | 0 | 1 | ? | ? | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | ? | ? | 1 | 0 | 0 | 0 | 0 | 3 | ? | 0 |
| H. lehmanniana | * | ? | 0 | 0 | 0 | 0 | 1 | 1 | ? | ? | ? | 1 | 0 | 0 | ? | 0 | ? | 1 | 0 | ? | 1 | 0 | ? | 0 | 0 | 1 | ? | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | ? | 0 | 0 |
| H. linearis | 0 | ? | 0 | 0 | 0 | 0 | ? | ? | ? | 0 | 0 | 0 | 1 | 0 | ? | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | ? | 0 | 1 | 0 | 0 | 0 | 0 | ? | 1 | 1 | 1 | 0 |
| H. lorea | 1 | ? | 0 | 0 | 1 | * | 0 | 1 | ? | ? | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 2 | 0 | 2 | 1 | ? | 2 | 0 | ? | 0 | ? | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | ? | 0 | 0 | 1 | 1 | 1 | 1 | 0 |
| H. megadenia | ? | ? | 0 | 0 | * | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | ? | 0 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 |
| H. megalosperma | 1 | 1 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 3 | 1 | 0 |
| H. multilineata | 0 | ? | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | ? | 2 | 1 | 2 | 1 | 0 | 2 | 0 | 0 | 1 | ? | ? | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | ? | ? | 1 | 0 | 0 | 0 | 0 | 0 | ? | 0 |
| H. nitida | 1 | 1 | 0 | 0 | 0 | 0 | ? | ? | 1 | ? | 0 | 0 | 1 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 2 | ? | 0 |
| H. obliqua | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 2 | 0 | 0 | 0 | 0 | 1 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 3 | ? | 1 |
| H. orthorrhyncha | 1 | ? | 0 | 0 | * | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 2 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | ? | 0 |
| H. pandanicarpa | 0 | 0 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 3 | ? | 1 |
| H. persiehana | 1 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 2 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 |
| H. petiolaris | * | 1 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | ? | 0 | 1 | 2 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 3 | 0 | 0 |
| H. platysperma | ? | 0 | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 2 | 0 | 1 | 1 | 0 | 1 | 2 | 2 | 0 | 1 | 0 | ? | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 1 | 3 | ? | 1 |
| H. propinqua | 0 | ? | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | ? | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | ? | 0 | 0 | 1 | 0 | 0 | 0 | 1 | ? | 0 |
| H. prostrata | 1 | 0 | 0 | 1 | 0 | 0 | ? | ? | 1 | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | ? | 0 | 2 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 1 | ? | 0 |
| H. purpurea | 1 | 1 | 0 | 0 | 1 | * | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 1 | 2 | ? | 0 | 1 | 0 | ? | 2 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 3 | 1 | 0 |
| H. pycnoneura | 0 | ? | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | ? | 0 | ? | 2 | 1 | 1 | 1 | 0 | 2 | 0 | 0 | 1 | ? | ? | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | ? | ? | 1 | 0 | 0 | 0 | 0 | 3 | ? | 0 |
| H. recurva | 0 | ? | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | ? | 0 | 2 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 1 | 1 | ? | ? | 0 |
| H. ruscifolia | 1 | 1 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | 0 | 2 | 0 | 0 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 2 | 1 | 0 |
| H. salicifolia | 0 | 0 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 0 | ? | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 |
| H. sericea | 0 | ? | 0 | 0 | 1 | 0 | 0 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 2 | ? | 0 |
| H. strumosa | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 3 | 1 | 1 |
| H. subsulcata | 0 | ? | 0 | 0 | 1 | 0 | ? | ? | ? | ? | ? | 1 | 0 | 0 | ? | 0 | ? | 2 | 1 | ? | 1 | 0 | ? | 0 | 0 | 1 | ? | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | ? | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 |
| H. teretifolia | * | ? | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 |
| H. trifurcata | 0 | 0 | 0 | 0 | * | * | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 2 | 1 | 1 | 0 | 0 | 0 | 2 | 2 | 1 | 1 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | 0 | 0 |
| H. trineura | 1 | 0 | 0 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | ? | ? | 0 | 0 | 2 | 0 | 2 | 1 | 0 | 2 | 0 | 0 | 1 | ? | 1 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | ? | 0 | 0 | 0 | 1 | ? | 0 |
| H. ulicina | 0 | ? | 1 | 0 | 0 | 0 | ? | ? | 0 | 0 | 1 | 1 | 0 | 0 | ? | 0 | ? | ? | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | ? | 0 | 1 | 0 | 0 | 0 | 0 | 2 | 0 | 0 |
| H. verrucosa | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | ? | ? | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 2 | ? | 0 | 1 | 0 | ? | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 |
| H. victoria | 0 | ? | 0 | 1 | 0 | 0 | ? | ? | 1 | ? | 1 | 1 | 0 | 0 | 0 | 0 | ? | 1 | 0 | 1 | 1 | 1 | 2 | 0 | 0 | 1 | 0 | 1 | 1 | ? | 0 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 2 | ? | 0 |
| Grevillea juncifolia | ? | ? | 0 | 0 | 0 | * | ? | ? | 0 | 0 | 0 | ? | 1 | 1 | 0 | 1 | 0 | 2 | 0 | 2 | 1 | 1 | 2 | ? | ? | ? | ? | 1 | 0 | 0 | 1 | 0 | ? | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | 0 | ? | 3 | ? | 0 |
Inferring phylogeny
We used the Akaike information criterion (Akaike, 1974) in MrModelTest 1.1b (available from J. A. A. Nylander, Uppsala University, Sweden, http://www.abc.se/∼nylander/) to select an adequately parameter-rich model of nucleotide substitution for each of the DNA regions. For the morphological data, we used the standard discrete model of MrBayes 3.1.2, which is based on the ideas of 39. These models then were used for their respective partitions in the Bayesian analyses in MrBayes (63). We unlinked the sampling of state frequencies, substitution rates, the gamma shape parameter, and the proportion of invariant sites for each DNA partition, and we linked branch lengths. We had MrBayes spawn two Markov Chain Monte Carlo (MCMC) runs when each partition was analyzed separately, and each of these MCMC runs ran four chains for 10 × 106 generations. The program spawned three MCMC runs that each ran four chains for 20 × 106 generations when the two partitions (the introns and the exons) of a waxy locus were analyzed together and four MCMC runs that each ran 12 chains for 40 × 106 generations when more than two partitions were analyzed together. MrBayes sampled every 1000th generation in all analyses. All MrBayes analyses were run in parallel on the Florida State University High Performance Computing Cluster. These numbers of MCMC chains and generations resulted in potential scale-reduction factors (25) for all parameters within at least two units of 1 (i.e., <3; and typically within 0.1 of 1). We constructed a majority-rule consensus of the trees sampled after the burn-in in PAUP* v4.0b10 (69) to produce a tree where the posterior probabilities were ≥95% (or ≥50%, when morphological characters were analyzed alone) for internal nodes.
For comparison, we also ran a maximum-parsimony search for the morphological data in PAUP* (69). PAUP* performed a full heuristic search, tree-bisection-reconnection branch swapping, and 10 random-addition replicates. Because this resulted in more trees than could be stored in memory, we saved 100000 of the most parsimonious trees, created a strict consensus, and ran another search, using the same heuristic search settings, looking for equally parsimonious trees that were incongruent with this strict consensus. The morphological data also were bootstrapped in 1000 replicates by the same search strategy as for the maximum-parsimony search but with a maximum number of trees set at 10000.
Inferring ages
We used Thorne's Bayesian method (72; 71) to infer the age of divergence events using the six coding regions because the taxa outside subtribe Hakeinae were not represented in the intron partitions. Thorne's approach relaxes the assumption of a constant rate of nucleotide substitution over time by modeling rate change with temporal autocorrelation (26). We explored the robustness of our results to alternative prior distributions for the rate-change parameter ν (described by brownmean and brownsd) and the rate at the root node (described by rtrate and rtratesd), as described later. 65 demonstrated a strong correlation between the ages of divergence events inferred by Thorne's method and penalized likelihood (64) and the uncorrelated log-normal method implemented in BEAST (19), respectively, in a family level study of the Proteaceae.
We used 81 baseml program in PAML version 3.13 and Thorne's estbranches and multidivtime programs in the way described by 50. The tree topology inferred from the combined molecular and morphological data (the “total-evidence tree”; 36) was used for the final step with multidivtime. We used a mean of 1.15 (representing 1.15 × 108 yr) for the prior distribution of time separating the most recent common ancestor (MRCA) of the ingroup (for this purpose, that of Platanaceae and Proteaceae) from the present (rttm). It is the mean of two age estimates for that node in a broad-scale dating study of the eudicots by 2. The standard deviation of this prior distribution (rttsd) was set to 0.15, a span that includes the earliest fossils that can be assigned to the stem of the Platanaceae, as described later. We used a mean of 0.0425 for the prior distribution of the rate of molecular evolution at the ingroup root node (rtrate). This value is one-half of the median distance between Plantanus L. and each taxon in the sister group (0.0978), as determined with the uncorrected-p parameter in PAUP*4.0b10 (69), divided by our prior estimate of time between the root of the ingroup and the present (1.15). The standard deviation of this prior distribution (rtratesd) was set to 0.0425 (equal to the mean; see 78, for a justification of making the standard deviation of this prior large relative to the mean). The “bigtime” parameter, which specifies what we consider to be a maximum age for the ingroup (in this case, a group composed of members of the Proteaceae and Platanaceae), was set at 1.3, equivalent to 130 Ma, older than the first appearance of the eudicots in the fossil record (45). The mean (brownmean) and standard deviation (brownsd) for the Brownian motion constant ν were both set to 1.0, as was the prior for the times of the interior nodes given the time of the root, according to recommendations made by Thorne in the multidivtime readme file. We left other parameters dealing with proposals in the MCMC at their default values (newk = 0.1, othk = 0.5, thek = 0.5).
Three constraints were used in the analyses. The maximum age for the MRCA of Proteaceae and Platanaceae was set at 125 Ma, the approximate time of appearance in the fossil record of the tricolpate pollen of the eudicots (45), of which Proteales is an early diverging branch. It has been used elsewhere as a calibration for the MRCAs of the eudicots (e.g., by 7), the Proteales and Sabiaceae (65), and the Proteaceae and Platanaceae (50). The minimum age for the MRCA of the Proteaceae and Platanaceae was set at 99.6 Ma, the upper boundary of the Albian (27), the earliest period from which fossil inflorescences have been recovered that can be assigned to Platanaceae (see, e.g., 10) on the basis of shared derived features (reviewed by 11; 2). It is slightly older than the occurrence of the oldest fossil that can be assigned to the Proteaceae from its position in parsimony analyses (65): Triorites africaensis (74; 16) from the Upper Cenomanian to Turonian (ca. 94 Ma) of Senegal and Gabon. The minimum age for the MRCA of Embothrium J.R.Forst. & G.Forst. and its sister, Telopea R.Br., was set to 35.4 Ma on the basis of the occurrence of Granodiporites nebulosus at that time (41; 65). Granodiporites nebulosus was resolved as sister to Embothrium by 65. 4 and 65 both used this fossil calibration in their studies of the family. 4 used a second calibration within tribe Embothrieae (the fossil pollen species Propylipollis ambiguus), but we do not use it here because 65 demonstrated that it cannot be assigned such precision with confidence. 65 also demonstrated that Hakeidites martinii, a pollen fossil previously assigned to Hakea by 35, cannot be assigned this precision with confidence.
We explored the sensitivity of the inference of ages to assumptions made in the dating analysis by systematically varying the combinations of the rtrate values of 0.00425, 0.0425 (the original value), and 0.425 and the brownmean values of 0.1, 1 (the original value), and 10; the standard deviation of each prior was set equal to the mean. A similar set of sensitivity tests was performed by 78, 7, and 50. We also determined the effect of removing the minimum age constraint for the MRCA of Embothrium and Telopea in one analysis and, in another, the minimum and maximum age constraints for the MRCA of Proteaceae and Platanaceae. In a final test of sensitivity, we increased the bigtime parameter to 2.6 and 13, in turn, while using the original values for rtrate and brownmean.
Coding pollen vectors
30 examined the distribution of SNDs from 25 species of Hakea for which floral visitors have been observed or postulated. They found all known or putative insect-pollinated species to have SND <13 mm and all known or putative bird-pollinated species to have SND >13 mm. They presented evidence of why SND is predictive by noting that the shortest bill length of the principal bird pollinators in Australia (the honeyeaters, family Meliphagidae) is 12 mm (58) and that insect pollinators in Western Australia rarely exceed a body length of 15 mm. We took SND from 30 for species included in that study; approximated SND for most of the remaining species from straightened pistil lengths given by 5, taking into account whether the pistil is curved; and measured SND from herbarium specimens when a species’ pistil was slightly longer than 13 mm and known to be curved. Although both the bill length of the honeyeaters and body length of insect pollinators might have differed considerably in the past from today's size ranges, 30 demonstration of a significant evolutionary correlation between inflorescence color (red is often, though not exclusively, associated with bird pollination) and SND size class in Hakea suggests that SND is a good predictor of pollinator class throughout the history of the genus.
Inferring ancestral states
We used two methods to determine ancestral states on the phylogeny. The Ancestral State Reconstruction Packages for Mesquite (43) in Mesquite (44) inferred character state changes from 22 parsimony on the chronogram inferred on the basis of our preferred set of prior assumptions. In the other method, MultiState in BayesTraits (57) calculated the maximum likelihood of the data when the MRCA of Hakea in the same chronogram was held as one state, then the other, using the “fossil” command. We used a log-likelihood difference of at least 2 as our threshold of significance, which is a rule of thumb used in the same way elsewhere (e.g., 56), following 20. We used a two-parameter model in which the evolutionary rates between the two states could be unequal because the log-likelihood of the data was significantly better (P = 0.0144) when we used that model (–28.526399) compared with the one-parameter model of equal rates (–31.519877), as determined using a likelihood ratio test with 1 degree of freedom. We excluded the outgroup taxa in these analyses to avoid introducing biases in estimates of evolutionary rate that might occur when a single taxon is used to represent a much larger group (as discussed, e.g., by 52). BayesTraits cannot calculate likelihood values for trees with polytomies, so we had Mesquite randomly resolve the two polytomies in the chromatogram and give the new branches a length of 0.00001 (the minimum length accepted by BayesTraits; the equivalent of 1000 yr).
Tests of evolutionary correlation, ordered evolution, and trait contingency
Because our results suggested that the morphology-based phylogeny of Hakea by 6 should be rerooted, we took the additional step to reanalyze the data on pollinators and floral cyanide from 30 following the same approach (from 56) and protocol of analysis that they did using the BayesDiscrete part of BayesTraits (57), with the original rooting and the new rooting. Specifically, we used the data and alternative rootings to determine whether evidence exists for a correlation between SND and floral cyanide concentrations, for the contingency of the evolution of bird pollination on high floral cyanide levels, and for greater rates from insect pollination + low cyanide concentrations to insect pollination + high cyanide concentrations than from the former to bird pollination + low cyanide concentration. Because their Table S1 did not include H. clavata Labill. or H. megalosperma Meisn., we excluded these from our analyses. 30 considered mean values for their floral cyanide assay ≤5 Feigl-Anger units (FA) to represent low cyanide concentrations and mean FA values ≥7.5 to represent high concentrations. The FA values for three species, H. erecta Lamont (mean = 5.5), H. recurva Meisn. (5.4), and H. smilacifolia Meisn. (6.8), fall between these thresholds. For our analyses, we considered the first two to have low floral cyanide concentrations and the last to have high concentrations. We performed the analyses with branch lengths equal (all set to 1), as done by 30. For the analysis, we used the topology shown in Figure 3 of 30. 30 referred to analyses performed on 10 trees that represented the most divergent possible resolutions of polytomies in the 6 topology. We did not attempt to recreate the other nine trees as well because our interest was not in exhaustively reanalyzing 30 results but rather in whether a rerooting of 6 topology made a difference.
RESULTS
Defining substitution characters
The sampled DNA regions are described in Table 3. We generated up to 9322 aligned positions for 75 taxa—ca. 700 kilobases of data in total. Each of the 75 species is represented in each set of coding DNA (which excludes the rpl16 and waxy introns) with seven exceptions (Appendix 1). Approximately 75% of the sites in the waxy intron alignments could be unambiguously aligned across subtribe Hakeinae; the remaining intron sites were not used in our study. Alignments of the genes matK, atpB, and ndhF represent the complete sequences for those genes without bordering regions. The alignment of the rpl16 intron does not include 36 nucleotides from the 5′ end and seven nucleotides from the 3′ end of the intron. Alignment of waxy exons 7–10 includes only 24 nucleotides on the 3′ end of exon 7 and only 171 nucleotides on the 5′ end of exon 10; exons 8 and 9 are complete. The DNA alignments are available as Appendix S1 (see Supplemental Data with the online version of this article).
| Region | Aligned length | Parsimony-informative nucleotide positions (% of aligned positions within Hakea) | Greatest sequence distance within Hakea | Substitution model selected on the basis of the Akaike information criterion |
|---|---|---|---|---|
| matK gene | 1566 | 64 (4.1%) | 0.01971 | GTR+I+G |
| atpB gene | 1440 | 32 (2.2%) | 0.01406 | GTR+I+G |
| ndhF gene | 2280 | 113 (5.0%) | 0.02369 | GTR+I+G |
| rpl16 intron | 1003 | 44 (4.4%) | 0.02643 | GTR+I+G |
| waxy 1 exons | 616 | 42 (6.8%) | 0.04181 | GTR+I+G |
| waxy 1 introns | 321 | 77 (24.0%) | 0.07220 | HKY+G |
| waxy 2 exons | 616 | 31 (5.0%) | 0.03496 | GTR+I+G |
| waxy 2 introns | 277 | 52 (18.8%) | 0.09438 | HKY+G |
| PHYA gene | 1169 | 60 (5.2%) | 0.03219 | HKY+I+G |
| Total | 9288 | 515 (5.5%) |
Coding morphological characters
Table 1 shows the morphological characters sampled, and table 2 shows the character-state data for each taxon. Including the three additional species of Grevillea used by 6; G. anethifolia R.Br., G. glauca Banks & Sol. ex Knight, and G. hakeoides Meisn.) in the morphological analysis did little to change the results, and because they were not sampled in the DNA data, they are not included in the analyses presented here.
Inferring phylogeny
Table 3 shows the models chosen as adequately parameter-rich for each set of data. The trees inferred from the individual cpDNA regions showed no topological incongruities (Appendix S2). Fig. 1 presents the trees for the four cpDNA regions combined and the individual nDNA regions; Fig. 2 presents the results of the parsimony and Bayesian analyses of the morphological data. Few incongruities exist between each of these topologies and the total-evidence topology (Fig. 3): the combined cpDNA topology is incongruous in four places; the waxy locus 1 topology in four; the waxy locus 2 topology in six; the PHYA topology in eight; the parsimony result from the morphological data in one; and the Bayesian result from the morphological data in 10.

Phylogenies inferred for subtribe Hakeinae when the four regions from DNA of the chloroplast and other plastids (cpDNA) were analyzed together and each of the nuclear DNA regions were analyzed separately. Only branches with 95% posterior probability or greater are shown. The asterisk (*) marks branches incongruent with the phylogeny inferred when all seven DNA regions and the morphological data were analyzed together (Fig. 3; the total-evidence phylogeny). The waxy gene trees were inferred from both exons and introns. The circumscription of Clade A and Clade B of the total-evidence phylogeny are indicated with light- and dark-gray boxes and labels. Because Clade A and Clade B are not monophyletic in all trees, they are labeled with quotation marks. Although 16 additional outgroup taxa were included in these analyses, we focused only on the results for subtribe Hakeinae to save space. The vertical dotted line at the base of the waxy locus 1 phylogeny indicates that lineages outside of subtribe Hakeinae appeared in the polytomy that includes Buckinghamia and Opisthiolepis.

Phylogenies inferred for the morphological data by parsimony (MP) and Bayesian (Bayes) methods, as indicated. The MP tree is the strict consensus of the most parsimonious trees for the data; bootstrap frequencies (as percentages) greater than 50% are given for that tree. The posterior probability of each branch inferred by the Bayesian method is given. The asterisk (*) marks branches incongruent with the phylogeny inferred when all seven DNA regions and the morphological data were analyzed together (Fig. 3; the total-evidence phylogeny). The circumscription of Clade A and Clade B of the total-evidence phylogeny are indicated with light- and dark-gray boxes and labels. Because Clade A and Clade B are not monophyletic in these trees, they are labeled with quotation marks.

Phylogeny inferred for subtribe Hakeinae from all the data with the membership of each species in 5 informal groups. The 12 informal groups that are not resolved as monophyletic in the phylogeny are in bold. The number of species sampled by this study and the total number of species in each group are given in parentheses. Although 16 additional outgroup taxa were included in this analysis, we focused only on the results for subtribe Hakeinae to save space. Branches with a “less-than” symbol (<) have a posterior probability between 50 and 95%; all other branches have a posterior probability ≥95%. Clade A and Clade B are informal names used in the text. The Multilineata clade recognized by 6 and shown here to be monophyletic is indicated.
The great number of incongruities with the PHYA topology arises because of an unexpected clade composed of members from both clades formed at the basal split of the total-evidence tree (Clades A and B; Fig. 3). The other two clades that share a polytomy with this clade in the PHYA topology are otherwise identical to Clades A and B, and they are labeled as such using quotation marks. Three processes are generally recognized as leading to incongruity among gene (or “linkage partition,” such as is the case for cpDNA) trees (53; 17; 42; 68): (i) lateral gene transfer, including hybridization and introgression; (ii) gene duplication (which can lead to paralogous comparisons); and (iii) deep coalescence (i.e., lineage sorting). We were unable to determine the cause of this incongruity: hybridization is unknown among extant populations of the six species in the anomalous clade; we did not observe two bands in our gel electrophoreses of PHYA for any species of Hakea; we did not find a second locus in the 10 colonies that we sequenced from the cloning of each of the six species involved; and the branch lengths between H. lorea and H. multilineata are not particularly long, as might be expected if they represent a period of time that includes the basal divergence in Hakea. Because the relationships inferred in the combined DNA and total-evidence trees did not depend on inclusion of this anomalous clade (but branch lengths did), we excluded the PHYA data for these six species in the combined DNA and total-evidence analyses (following 38).
The basal split in the morphology-based results is between H. lorea and the remaining sampled species. Although this rooting is congruent with the results of 6, the other early branching events are not. Few branches in the morphology-based results were well supported though (e.g., there was just one branch with bootstrap support >70% and two with posterior probabilities >90%), and all but one of the 10 incongruities between the Bayesian result and the total-evidence phylogeny had <90% posterior probabilities in the former. 6 did not determine support values for branches, so we cannot compare the results of the two studies in this regard.
The total-evidence topology (Fig. 3) is identical to that found when the seven DNA regions were analyzed in combination without the morphological data. Of the 19 informal groups of 5 from which two or more species were sampled in our study, 12 are resolved as nonmonophyletic and seven as monophyletic in the total-evidence tree (Fig. 3). Of the five informal groups of 5 used by 30 to identify closely related species for comparison (Hanley et al. table 2; their species group “Sulcata” is not separate from the Ulicina group in 5, and so is not counted here), four of them are not monophyletic in our total-evidence tree. Five of the six multigroup clades named by 6 are not monophyletic on the phylogeny from our analysis of molecular and morphological data combined, but the largest of these, the Multilineata clade, is monophyletic and makes up most of Clade B (Fig. 3).
Inferring ages
We inferred the age of the crown group of Hakea to be 9.6 Ma (= mean; 95% credible interval = 6.4–14.0) using our preferred set of parameter values and all three age constraints. The mean inferred age varied from 8.1 Ma to 17.4 Ma as we varied values of rtrate, rtratesd, brownmean, and brownsd. The highest age estimates were inferred when the prior of the Brownian motion constant ν was set to 0.1; this is the most “clock-like” value assumed for that variable, an assumption that is perhaps inappropriate given the distribution of branch lengths (see phylogram, Fig. 4c). The mean of all other age estimates for the crown group was <10 Ma. For example, removing the minimum age constraint (35.4 Ma) for the MRCA of Embothrium and Telopea resulted in a mean inferred age of 7.7 Ma (4.8–11.3) for the crown group of Hakea, and removing the minimum and maximum age constraint for the MRCA of Proteaceae and Platanaceae resulted in a mean inferred age of 9.6 Ma (6.3–13.8). Increasing the bigtime parameter did not substantially change the ages inferred by the method. The second replicate for each of these analyses resulted in mean and boundaries of the 95% credible interval values within 0.3 Ma of those of the first replicate. Using our preferred set of parameter values and all three age constraints, we inferred the age of the stem group of Hakea to be 12.8 Ma (8.8–18.0).

Chronogram and phylogram for order Proteales with ancestral state reconstructions for pollination in Hakea. (A) Chronogram. The position of each node on the x-axis of the chronogram is at the mean inferred age for that node with three calibration points and our preferred set of prior assumptions (mean and standard deviation of the prior distribution of the rate of molecular evolution at the ingroup root node, rtrate and rtratesd = 0.0425 and mean and standard deviation of the prior distribution of the Brownian motion constant, brownmean and brownsd = 1). 95% credible intervals are given for the crown group of Hakea and the most recent common ancestor of Hakea and the sampled species of Grevillea and Finschia. (B) Phylogram. The phylogram shows the mean branch lengths from the posterior probability density for the analysis with the six coding DNA regions. (C) Ancestral reconstruction of stigma–nectary distance, an indication of insect (<13 mm) and bird pollination (>13 mm). Branch shading indicates the results of Fitch parsimony reconstruction. Branches on which multiple states are inferred are shaded gray. The likelihood of the data is significantly better when the MRCA of Hakea is held as bird pollinated than when it is held as insect pollinated.
Coding pollen vectors
Of 19 species included in both our study and that of 30, nine have SND <13 mm and 10 SND >13 mm. Of the 36 species for which we needed to estimate SND from pistil lengths from 5, 27 had pistil length <13 mm and nine had pistil length >13 mm. Of those with pistil length >13 mm, two included 13 mm in the bottoms of their ranges: H. pandanicarpa R.Br. (13–16 mm) and H. grammatophylla F.Muell. (13–16 mm). The flowers of H. pandanicarpa have a strongly incurved pistil. We measured a mean SND of 11.9 mm (range 9–11.5 mm) for H. pandanicarpa from a single flower from each of five specimens at the National Herbarium of New South Wales (NSW). We consider this result and the strong, sweet scent of its flowers (82) to be evidence that it is insect pollinated. The flowers of H. grammatophylla have a bent, but not strongly incurved, pistil. We measured a mean SND of 14.9 mm (range 13–17) for H. grammatophylla from a single flower from each of five specimens at NSW. We consider this and its deep pink flowers in large, exposed inflorescences to be evidence that it is bird pollinated. These codings of pollinator for the species are consistent with those proposed by 6 in cases of overlap with that study, with the exception that 6 considered the Corymbosa group (represented by H. victoria and H. corymbosa here) to be mammal, rather than bird, pollinated (contra 30).
Inferring ancestral states
Parsimony inference of ancestral SNDs resolve the MRCA of Hakea to have been bird pollinated (to have had SND >13 mm) with either nine shifts to insect pollination (SND <13 mm) or seven shifts to insect pollination and two back to bird pollination (Fig. 4b). BayesTraits calculated the rate of change from bird to insect pollination is 10.752508, and the rate of change from insect to bird pollination is 0. The log-likelihood of the data is significantly different when the MRCA of Hakea is held as bird pollinated (–28.526399) as opposed to insect pollinated (–31.447590).
Tests of evolutionary correlation, ordered evolution, and trait contingency
Using 30 original data and tree (the topology of their Fig. 3), we found both the correlation between SND and floral cyanide concentrations (χ2 = 12.92, df = 4, P = 0.012) and evidence that the evolution of bird pollination (SND >13 mm) is contingent on high cyanide levels (contingent evolution; χ2 = 5.04, df = 1, P = 0.025) to be significant, but the evidence for ordered evolution (a greater evolutionary rate from insect pollination + low cyanide concentrations to insect pollination + high cyanide concentrations than to bird pollination + low cyanide concentrations; χ2 = 0.000002, df = 1, P = 0.999) not to be. Oddly, the P value we found for the test of ordered evolution (P = 0.999) was not within the range that 30 report for the test with each of the 10 tree topologies (P = 0.03–0.20), but our results for the correlated evolution and contingent evolution tests are similar to theirs. After rerooting their topology to make the ingroup rooting consistent with what we found, evidence for correlated evolution between the two characters remained significant (χ2 = 12.54, df = 4, P = 0.014), and evidence of ordered evolution remained nonsignificant (χ2 = 0.000002, df = 1, P = 0.999), but the evidence for contingent evolution became less significant (χ2 = 2.31, df = 1, P = 0.129).
DISCUSSION
Our study provides a well-supported backbone phylogeny for Hakea based on both morphological data and, for the first time, molecular data, and it demonstrates the critical role that assumptions about the phylogeny and age of a lineage can play in evolutionary studies. We show that the inferred phylogeny of Hakea changed substantially when molecular data were added to the morphological data of 6 and that the result profoundly changed the strength of evidence for a central conclusion of 30 study of the evolution of pollination systems in the genus: that bird-pollinated taxa had evolved multiple times from insect-pollinated ancestors.
Phylogenetic congruence with previous taxon circumscriptions
We inferred that the basal split in Hakea produced two clades (Fig. 3), each with moderate morphological cohesiveness. We estimate that the larger clade (Clade A) contains 96 species, on the basis of memberships in the informal groups of 5. All sampled members of each informal group, whether it is monophyletic or not, appear in either Clade A or Clade B but not both. Clade A is characterized by obscure leaf venation and discoid pollen presenters, with the exception of the Lissocarpha and Varia groups, which have conical pollen presenters. We estimate that Clade B contains about half as many species (53) as Clade A. It is composed primarily of the Multilineata clade recognized by 6; Fig. 3) and is characterized by prominent leaf venation and conical pollen presenters. Determination of the derived state in each of these two characters will have to await clarification regarding the sister group to Hakea, since it might not be the whole of Grevillea and Finschia. The two clades formed at the basal split have not been recognized previously.
Other putative multigroup clades of 6 and 12 of the 19 informal groups of 5 sampled multiple times in our study are nonmonophyletic in the total-evidence phylogeny (Fig. 3) and thus should not be considered monophyletic in future evolutionary studies. Although this result at first could be viewed as a striking incongruity with 5 study, we note that five of the 12 were not resolved as monophyletic in the cladogram of 6 either and that four other nonmonophyletic informal groups were represented by only one species in the 6 study, so their monophyly was not previously tested. Although branch-support values were not reported by 6, those we present (Fig. 2) suggest that the morphological data alone do not produce strong corroboration for the monophyly of the three remaining informal groups (relevant bootstrap values all <50%), so the incongruence is not worrisome.
The age of Hakea
The recent origin of the crown group of Hakea (ca. 10 Ma) is unexpected on the basis of both 30 expectations and comparisons with the ages of ecologically similar lineages. For example, in a study of diversification timing in four sclerophyll lineages in Australia, 12 found that much of the extant diversity of those groups arose from rapid diversification throughout an earlier period that corresponds mostly to the early and mid Miocene (25–10 Ma ago, according to 12). At that time, Australia's climate is thought to have become drier and more seasonal, as marked by a decrease, followed by disappearance, of Nothofagus pollen and a marked increase in pollen of Myrtaceae and other sclerophyll taxa in the fossil record (47). Although good phylogenetic evidence exists for the evolution of adaptations to fire as early as the Cretaceous–Palaeogene boundary (14), fire is thought to have become more regular in the late Miocene (47), a period between 11.6 and 5.3 Ma ago (27) that provides a reasonable context for the increase in follicle-wall thickness that is a synapomorphy for the genus Hakea. This conclusion is consistent with 6 view that the Miocene would have been a time conducive to the radiation of Hakea, had this lineage originated by then.
Our Bayesian age estimates were robust against large variations in our prior assumptions, and our age estimate for the stem group of Hakea (mean = 12.8 Ma, 95% credible interval = 8.8–18.0 Ma; Fig. 4a) is similar to that found in the family-wide study by 65; mean = 14.1 Ma, 95% credible interval = 8.3–21.52 Ma) when they also used Thorne's Bayesian method (72; 71). Using a method in which autocorrelation in rates in adjacent branches is not assumed (19), they estimated the stem age to be only slightly older (mean = 16.2 Ma; 95% credible interval = 9.2–23.3 Ma). 65 use of a single OTU for Hakea means that they did not infer the age of its crown group. Another recent, family-wide dating study (4) did not include a representative of Hakea. If Grevillea proves to be paraphyletic and we failed to sample the sister lineage of Hakea in the current study, then our age estimate for the stem group of Hakea will prove to be higher, not lower, than justified.
The evolution of pollination
On our phylogeny of Hakea, parsimony and likelihood approaches reconstructed the MRCA of the genus to have had an SND that occurs in extant species that use birds as pollinators (>13 mm), rather than one that corresponds to insect pollination (<13 mm) as found in 30 more restricted geographic and taxonomic sampling. The principal bird pollinators of Hakea, the honeyeaters (6; 30), appear to have been present in Australia from at least the Miocene (8), making reasonable the inference of bird pollination in a 10-Ma-old MRCA of Hakea.
Our parsimony inference of multiple shifts from bird to insect pollination, rather than the reverse, and our maximum likelihood inference of a rate of 0 from insect to bird pollination (as opposed to 10.8 for the reverse) make Hakea only the third documented example of this kind of transition that we can find in the literature on pollinator shifts. However, that literature is heavily biased toward New World examples, which have focused on taxa in which the only bird-pollinated plants are those associated with hummingbirds. Ours is only the third phylogenetic study that has examined transitions between insect and honeyeater pollinators, all of which have raised the possibility of shifts from bird to insect pollination. More phylogenetic results like ours are required before such a pattern can reasonably be said to demand explanation, but striking differences in the foraging behavior of these two bird clades suggest that such a pattern might be expected. Hummingbirds, which hover when feeding, can be expected to show a finer degree of precision in their interactions with their food plants than honeyeaters, which perch to feed. Moreover, honeyeater-pollinated plants tend to be visited by a wide range of bird species of dramatically variable size, which themselves feed from a wide range of plant species (23). Even species of Hakea from Australia's arid center are exposed to a fauna of 11 meliphagid species, ranging in length from 10 to 28 cm (67). This lack of specialization compared with hummingbird-pollinated plants implies the possibility of greater flexibility in shifting to generalist pollination and thence to more specialized insect pollination than for hummingbird-pollinated taxa.
Our inability to replicate 30 “ordered evolution” and, with a rerooting of the ingroup, “contingent evolution” results suggest that 30 conclusion that bird pollination only arose where insect-pollinated ancestors were already producing high concentrations of floral cyanide is incorrect. That said, new attention might be productively shifted to groups diverging earlier in the phylogeny than the crown group of Hakea because the bird pollination primitive to the crown group can be assumed to have arisen at some time leading up to it, and neither our study nor that of 30 examined that pollination shift. 23 estimated that 150 of the 357 species of Grevillea are bird pollinated.
The divergence events leading up to the MRCA of Hakea and the species of Grevillea sampled in our study might reasonably be expected to have produced other lineages of Grevillea, given that our results and the monophyly of Grevillea would otherwise imply that the diversification that produced the third largest genus in the Australian flora (357 species) occurred over a period of just 13 Ma. A young age for Grevillea would be at odds with the results of 40, who, in a review of published studies of clade ages, concluded that the Australian flora appears to consist largely of slowly diversifying, ancient lineages (>15 Ma old). On the other hand, it would support the hypothesis of 32, 33) that a period of slower speciation was followed, in the late Tertiary and Quaternary, by an acceleration of speciation in southwestern Australia (where large fractions of Hakea and Grevillea are endemic). Clearly, further sampling of Grevillea (Olde and Marriott 1995; Makinson 2000) is needed before hypotheses regarding the evolutionary interplay between pollination and florivory in subtribe Hakeinae (a now justifiably broader focus than just Hakea) and diversification rates in the subtribe can be adequately addressed.
Appendix 1
Taxa included in this study with their collection voucher (with the herbarium where it can be found in parentheses using the abbreviations of Index Herbariorum). Abbreviations: n/a = taxa for which the rpl16intron was not a target; (dash) = target regions that could not be sequenced.
| Chloroplast regions | Nuclear regions | |||||||
|---|---|---|---|---|---|---|---|---|
| Species | Voucher | matK gene | atpB gene | ndhF gene | rpl16 intron | PHYA gene | waxy locus 1 | waxy locus 2 |
| Fam. Nelumbonaceae | ||||||||
| Nelumbo lutea | Sarah Braun et al. s.n. (FSU) | EU642710 | EU642740 | EU642680 | n/a | EU642772 | EU649747 | EU649747 |
| Fam. Platanaceae | ||||||||
| Platanus occidentalis | Austin R. Mast s.n. (FSU) | EU642711 | EU642741 | EU642681 | n/a | EU642773 | EU649748 | EU649748 |
| Fam. Proteaceae | ||||||||
| Subfam. Bellendenoideae | ||||||||
| Bellendena montana | Greg Jordan s.n. (NSW) | JQ257173 | JQ257242 | JQ257311 | n/a | JQ257439 | — | — |
| Subfam. Persoonioideae | ||||||||
| Persoonia lanceolata | Peter H. Weston 2489 (NSW) | JQ257172 | JQ257241 | JQ257310 | n/a | JQ257438 | — | JQ313642 |
| Subfam. Symphionematoideae | ||||||||
| Symphionema montanum | Greg Jordan s.n. (NSW) | JQ257174 | JQ257243 | JQ257312 | n/a | JQ257440 | JQ313709 | JQ313643 |
| Subfam. Proteoideae | ||||||||
| Protea cordata | Austin R. Mast 475 (NBG) | JQ257175 | JQ257244 | JQ257313 | n/a | JQ257441 | — | — |
| Subfam. Grevilleoideae | ||||||||
| Tr. Roupaleae | ||||||||
| Roupala montana | Peter H. Weston 2038 (NSW) | EU642684 | EU642714 | EU642652 | n/a | EU642744 | EU649751 | EU649719 |
| Tr. Banksieae | ||||||||
| Banksia serrata | Austin R. Mast et al., 225 (WIS) | AY823169 | AY837794 | EU642656 | n/a | EU642748 | AY829483 | EU649723 |
| Tr. Macadamieae | ||||||||
| Macadamia integrifolia | FTG 355A (WIS) | AY823204 | AY837827 | EU642672 | n/a | EU642764 | EU649739 | AY829522 |
| Tr. Embothrieae | ||||||||
| Subtr. Lomatiinae | ||||||||
| Lomatia silaifolia | Peter H. Weston 2875 (NSW) | JQ257167 | JQ257236 | JQ257305 | n/a | JQ257433 | JQ313704 | JQ313637 |
| Subtr. Embothriinae | ||||||||
| Embothrium coccineum | NSW 863505D (NSW) | JQ257169 | JQ257238 | JQ257307 | n/a | JQ257435 | JQ313706 | JQ313639 |
| Oreocallis mucronata | Peter H. Weston 1927 (NSW) | JQ257171 | JQ257240 | JQ257309 | n/a | JQ257437 | JQ313708 | JQ313641 |
| Alloxylon flammeum | D. C. Godden 277 (NSW) | JQ257168 | JQ257237 | JQ257306 | n/a | JQ257434 | JQ313705 | JQ313638 |
| Telopea speciosissima | D. C. Godden 207 (NSW) | JQ257170 | JQ257239 | JQ257308 | n/a | JQ257436 | JQ313707 | JQ313640 |
| Subtr. Stenocarpinae | ||||||||
| Stenocarpus sinuatus | Peter H. Weston 2026 (NSW) | JQ257166 | JQ257235 | JQ257304 | n/a | JQ257432 | JQ313703 | JQ313636 |
| Strangea cynanchicarpa | Peter H. Weston 1975 (NSW) | JQ257165 | JQ257234 | JQ257303 | n/a | JQ257431 | JQ313702 | JQ313635 |
| Subtr. Hakeinae | ||||||||
| Buckinghamia celsissima | Carolyn L. Porter 419835 (NSW) | JQ257176 | JQ257245 | JQ257314 | JQ257373 | JQ257442 | JQ313710 | JQ313644 |
| Opisthiolepis heterophylla | Carolyn L. Porter 427629 (NSW) | JQ257177 | JQ257246 | JQ257315 | JQ257374 | JQ257443 | JQ313711 | JQ313645 |
| Grevillea juncifolia | Peter H. Weston 2913 (NSW) | JQ257178 | JQ257247 | JQ257316 | JQ257375 | JQ257444 | JQ313712 | JQ313646 |
| Finschia chloroxantha | Irian Jaya (NSW) | JQ257233 | JQ257302 | JQ257372 | JQ257430 | JQ257499 | JQ313767 | JQ313701 |
| H. arborescens | Paul Kennedy s.n. (NSW) | JQ257232 | JQ257301 | JQ257371 | JQ257429 | JQ257498 | JQ313766 | JQ313700 |
| H. archaeoides | Peter H. Weston 2803 (NSW) | JQ257207 | JQ257276 | JQ257346 | JQ257404 | — | JQ313741 | JQ313676 |
| H. auriculata | Austin R. Mast 613 (PERTH) | JQ257211 | JQ257280 | JQ257350 | JQ257408 | JQ257477 | JQ313745 | JQ313680 |
| H. baxteri | Peter M. Olde 04/278 (NSW) | JQ257217 | JQ257286 | JQ257356 | JQ257414 | JQ257483 | JQ313751 | JQ313686 |
| H. brachyptera | Austin R. Mast 617 (PERTH) | JQ257193 | JQ257262 | JQ257332 | JQ257390 | JQ257460 | JQ313727 | JQ313662 |
| H. bucculenta | Peter H. Weston 2802 (NSW) | JQ257206 | JQ257275 | JQ257345 | JQ257403 | JQ257473 | JQ313740 | JQ313675 |
| H. clavata | Peter H. Weston 2806 (NSW) | JQ257194 | JQ257263 | JQ257333 | JQ257391 | JQ257461 | JQ313728 | JQ313663 |
| H. commutata | Peter M. Olde 04/289 (NSW) | JQ257221 | JQ257290 | JQ257360 | JQ257418 | JQ257487 | JQ313755 | JQ313690 |
| H. conchifolia | Peter M. Olde 04/142 (NSW) | JQ257222 | JQ257291 | JQ257361 | JQ257419 | JQ257488 | JQ313756 | JQ313691 |
| H. constablei | Peter H. Weston 2857 (NSW) | JQ257210 | JQ257279 | JQ257349 | JQ257407 | JQ257476 | JQ313744 | JQ313679 |
| H. corymbosa | Peter H. Weston 2808 (NSW) | JQ257203 | JQ257272 | JQ257342 | JQ257400 | JQ257470 | JQ313737 | JQ313672 |
| H. cristata | Peter H. Weston 2816 (NSW) | JQ257186 | JQ257255 | JQ257325 | JQ257383 | JQ257453 | JQ313720 | JQ313655 |
| H. cucullata | Peter H. Weston 2807 (NSW) | JQ257202 | JQ257271 | JQ257341 | JQ257399 | JQ257469 | JQ313736 | JQ313671 |
| H. dactyloides | Peter H. Weston 2794 (NSW) | JQ257199 | JQ257268 | JQ257338 | JQ257396 | JQ257466 | JQ313733 | JQ313668 |
| H. divaricata | Peter H. Weston 2814 (NSW) | JQ257209 | JQ257278 | JQ257348 | JQ257406 | JQ257475 | JQ313743 | JQ313678 |
| H. drupacea | Peter H. Weston LCR 872062 (NSW) | JQ257195 | JQ257264 | JQ257334 | JQ257392 | JQ257462 | JQ313729 | JQ313664 |
| H. eriantha | Peter H. Weston 2797 (NSW) | JQ257187 | JQ257256 | JQ257326 | JQ257384 | JQ257454 | JQ313721 | JQ313656 |
| H. florida | Peter H. Weston 2810 (NSW) | JQ257196 | JQ257196 | JQ257335 | JQ257393 | JQ257463 | JQ313730 | JQ313665 |
| H. grammatophylla | Peter H. Weston 2825 (NSW) | JQ257213 | JQ257282 | JQ257352 | JQ257410 | JQ257479 | JQ313747 | JQ313682 |
| H. hastata | Peter M. Olde 04/211 (NSW) | JQ257223 | JQ257292 | JQ257362 | JQ257420 | JQ257489 | JQ313757 | JQ313692 |
| H. horrida | Peter M. Olde 04/199 (NSW) | JQ257219 | JQ257288 | JQ257358 | JQ257416 | JQ257485 | JQ313753 | JQ313688 |
| H. incrassata | Austin R. Mast 619 (PERTH) | JQ257188 | JQ257257 | JQ257327 | JQ257385 | JQ257455 | JQ313722 | JQ313657 |
| H. invaginata | Peter H. Weston 2809 (NSW) | JQ257227 | JQ257296 | JQ257366 | JQ257424 | JQ257493 | JQ313761 | JQ313696 |
| H. lasianthoides | Peter M. Olde 04/268 (NSW) | JQ257215 | JQ257284 | JQ257354 | JQ257412 | JQ257481 | JQ313749 | JQ313684 |
| H. laurina | Peter M. Olde 04/285 (NSW) | JQ257224 | JQ257293 | JQ257363 | JQ257421 | JQ257490 | JQ313758 | JQ313693 |
| H. lehmanniana | Peter M. Olde 04/216 (NSW) | JQ257229 | JQ257298 | JQ257368 | JQ257426 | JQ257495 | JQ313763 | JQ313698 |
| H. linearis | M. Pieroni 06/1 (PERTH) | JQ257228 | JQ257297 | JQ257367 | JQ257425 | JQ257494 | JQ313762 | JQ313697 |
| H. lorea | Peter H. Weston 2801 (NSW) | JQ257181 | JQ257250 | JQ257319 | JQ257378 | JQ257447 | JQ313715 | JQ313649 |
| H. megadenia | Greg Jordan s.n. (HO) | JQ257225 | JQ257294 | JQ257364 | JQ257422 | JQ257491 | JQ313759 | JQ313694 |
| H. megalosperma | Austin R. Mast 628 (PERTH) | JQ257179 | JQ257248 | JQ257317 | JQ257376 | JQ257445 | JQ313713 | JQ313647 |
| H. multilineata | Peter H. Weston 2805 (NSW) | JQ257201 | JQ257270 | JQ257340 | JQ257398 | JQ257468 | JQ313735 | JQ313670 |
| H. nitida | Peter M. Olde 04/279 (NSW) | JQ257220 | JQ257289 | JQ257359 | JQ257417 | JQ257486 | JQ313754 | JQ313689 |
| H. obliqua | Peter M. Olde 04/311 (NSW) | JQ257216 | JQ257285 | JQ257355 | JQ257413 | JQ257482 | JQ313750 | JQ313685 |
| H. orthorrhyncha | Peter H. Weston 2813 (NSW) | JQ257197 | JQ257266 | JQ257336 | JQ257394 | JQ257464 | JQ313731 | JQ313666 |
| H. pandanicarpa | Peter H. Weston 2817 (NSW) | JQ257198 | JQ257267 | JQ257337 | JQ257395 | JQ257465 | JQ313732 | JQ313667 |
| H. persiehana | Paul Kennedy s.n. (NSW) | JQ257231 | JQ257300 | JQ257370 | JQ257428 | JQ257497 | JQ313765 | JQ313699 |
| H. petiolaris | Peter H. Weston 2804 (NSW) | JQ257200 | JQ257269 | JQ257339 | JQ257397 | JQ257467 | JQ313734 | JQ313669 |
| H. platysperma | Peter M. Olde 04/126 (NSW) | JQ257218 | JQ257287 | JQ257357 | JQ257415 | JQ257484 | JQ313752 | JQ313687 |
| H. propinqua | Peter H. Weston 2792 (NSW) | JQ257183 | JQ257252 | JQ257322 | JQ257380 | JQ257450 | JQ313717 | JQ313652 |
| H. prostrata | Austin R. Mast 623 (PERTH) | JQ257185 | JQ257254 | JQ257324 | JQ257382 | JQ257452 | JQ313719 | JQ313654 |
| H. purpurea | Peter H. Weston 2799 (NSW) | JQ257212 | JQ257281 | JQ257351 | JQ257409 | JQ257478 | JQ313746 | JQ313681 |
| H. pycnoneura | Peter H. Weston 2826 (NSW) | JQ257208 | JQ257277 | JQ257347 | JQ257405 | JQ257474 | JQ313742 | JQ313677 |
| H. recurva | Austin R. Mast 620 (PERTH) | JQ257205 | JQ257274 | JQ257344 | JQ257402 | JQ257472 | JQ313739 | JQ313674 |
| H. ruscifolia | Peter H. Weston 2819 (NSW) | JQ257226 | JQ257295 | JQ257365 | JQ257423 | JQ257492 | JQ313760 | JQ313695 |
| H. salicifolia | Peter H. Weston 2822 (NSW) | JQ257189 | JQ257258 | JQ257328 | JQ257386 | JQ257456 | JQ313723 | JQ313658 |
| H. sericea | Peter H. Weston 2786 (NSW) | JQ257182 | JQ257251 | JQ257321 | JQ257379 | JQ257449 | JQ313716 | JQ313651 |
| H. stenophylla | Peter M. Olde 04/155 (NSW) | JQ257214 | JQ257283 | JQ257353 | JQ257411 | JQ257480 | JQ313748 | JQ313683 |
| H. strumosa | Peter H. Weston 2820 (NSW) | JQ257190 | JQ257259 | JQ257329 | JQ257387 | JQ257457 | JQ313724 | JQ313659 |
| H. subsulcata | Peter M. Olde 04/317 (NSW) | JQ257230 | JQ257299 | JQ257369 | JQ257427 | JQ257496 | JQ313764 | — |
| H. teretifolia | Peter H. Weston 2790 (NSW) | JQ257192 | JQ257261 | JQ257331 | JQ257389 | JQ257459 | JQ313726 | JQ313661 |
| H. trifurcata | Austin R. Mast 623.5 (PERTH) | JQ257191 | JQ257260 | JQ257330 | JQ257388 | JQ257458 | JQ313725 | JQ313660 |
| H. trineura | Peter H. Weston 2811 (NSW) | JQ257180 | JQ257249 | JQ257318 | JQ257377 | JQ257446 | JQ313714 | JQ313648 |
| H. ulicina | Peter H. Weston 2815 (NSW) | JQ257204 | JQ257273 | JQ257343 | JQ257401 | JQ257471 | JQ313738 | JQ313673 |
| H. verrucosa | Peter H. Weston 2812 (NSW) | JQ257184 | JQ257253 | JQ257323 | JQ257381 | JQ257451 | JQ313718 | JQ313653 |
| H. victoria | Austin R. Mast 076 (PERTH) | EU649719 | EU649719 | JQ257320 | EU649719 | JQ257448 | EU649719 | JQ313650 |